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Let X = (Xt)t>o be a transient diffusion process in (0, oo) with 
the diffusion coefficient a > and the scale function L such that 
X t — > oo as t — > oo, let It denote its running minimum for t > 0, and 
let 6 denote the time of its ultimate minimum I x . Setting c(i,x) = 
1 — 2L(x) / L(i) we show that the stopping time 

n=mi{t>0\Xt >/*(/«)} 

minimizes E(\0 — r\ — 8) over all stopping times r of X (with finite 
mean) where the optimal boundary /„ can be characterized as the 
minimal solution to 

m _ <r 2 (f(i))L'W)) f m c'S,y)[L{y)-L{i)] 

JU c{i,!{i))[L{}{{))-L{i)\j i a*(y)L'(y) V 

staying strictly above the curve h(i) = L _1 (L(i)/2) for i > 0. In par- 
ticular, when X is the radial part of three-dimensional Brownian 
motion, we find that 

r. = itf\t>o\ Xt ~ h > (p 



where ip = (1 + V / 5)/2 = 1.61 ... is the golden ratio. The derived re- 
sults are applied to problems of optimal trading in the presence of 
bubbles where we show that the golden ratio rule offers a rigorous 
optimality argument for the choice of the well-known golden retrace- 
ment in technical analysis of asset prices. 
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1. Introduction. The golden ratio has fascinated people of diverse inter- 
ests for at least 2400 years (see, e.g., [24]). In mathematics (and the arts) two 
quantities a and b are in the golden ratio if the ratio of the sum of the quan- 
tities a + b to the larger quantity a is equal to the ratio of the larger quantity 
a to the smaller quantity b. This amounts to setting (a + b)/a = a/b =: ip and 
solving ip 2 — ip — 1 = which yields ip = (1 + y/E)/2 = 1.61 . . . Apart from be- 
ing abundant in nature, and finding diverse applications ranging from archi- 
tecture to music, the golden ratio has also found more recent uses in technical 
analysis of asset prices (in strategies such as Fibonacci retracement repre- 
senting an ad-hoc method for determining support and resistance levels). 
Despite its universal presence and canonical role in diverse applied areas, 
we are not aware of any more exact connections between the golden ratio 
and stochastic processes (including any proofs of optimality in particular). 

One of the aims of the present paper is to disclose the appearance of the 
golden ratio in an optimal stopping strategy related to the radial part of 
three-dimensional Brownian motion. More specifically, denoting the radial 
part by X it is well known that X is transient in the sense that Xt — > oo 
as t — > oo. After starting at some x > 0, the ultimate minimum of X will 
therefore be attained at some time 9 that is not predictable through the 
sequential observation of X (in the sense that it is only revealed at the end 
of time). The question we are addressing is to determine a (predictable) 
stopping time of X that is as close as possible to 9. We answer this question 
by showing that the time at which the excursion of X away from the running 
minimum I and the running minimum I itself form the golden ratio is as 
close as possible to 9 in a normalized mean deviation sense. We consider this 
problem by embedding it into transient Bessel processes of dimension d > 2 
and in this context we derive similar optimal stopping rules. We also disclose 
further /deeper extensions of these results to transient diffusion processes. 
The relevance of these questions in financial applications is motivated by 
the problem of optimal trading in the presence of bubbles. In this context 
we show that the golden ratio rule offers a rigorous optimality argument 
for the choice of the well-known golden retracement in technical analysis of 
asset prices. To our knowledge this is the first time that such an argument 
has been found/given in the literature. 

The problem considered in the present paper belongs to the class of opti- 
mal prediction problems (within optimal stopping). Similar optimal predic- 
tion problems have been studied in recent years by many authors (see, e.g., 
[1, 2, 6-9, 14, 19, 25, 30, 39-41]). Once the "unknown" future is projected to 
the "known" present, we find that the resulting optimal stopping problem 
takes a novel integral form that has not been studied before. The appearance 
of the minimum process in this context makes the problem related to optimal 
stopping problems for the maximum process that were initially studied and 
solved in important special cases of diffusion processes in [11, 12] and [23]. 
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The general solution to problems of this kind was derived in the form of the 
maximality principle in [31]; see also Section 13 and Chapter V in [35] and 
the other references therein. More recent contributions and studies of related 
problems include [3, 15, 16, 20, 22, 26, 27, 29]. Close three-dimensional rel- 
atives of these problems also appear in the recent papers [10] and [43] where 
the problems were effectively solved by guessing and finding the optimal 
stopping boundary in a closed form. The maximality principle has been 
extended to three-dimensional problems in the recent paper [34]. 

Although the structure of the present problem is similar to some of these 
problems, it turns out that none of these results is applicable in the present 
setting. Governed by these particular features in this paper we show how the 
problem can be solved when (i) no closed-form solution for the candidate 
stopping boundary is available and (ii) the loss function takes an integral 
form where the integrand is a functional of both the process X and its 
running minimum /. This is done by extending the arguments associated 
with the maximality principle to the setting of the present problem and 
disclosing the general form of the solution that is valid in all particular cases. 
The key novel ingredient revealed in the solution is the replacement of the 
diagonal and its role in the maximality principle by a nonlinear curve in the 
two-dimensional state space of X and I. We believe that this methodology 
is of general interest and the arguments developed in the proof should be 
applicable in similar two/multi-dimensional integral settings. 

2. Optimal prediction problem. 1. We consider a nonnegative diffusion 
process X = (X t )t>o solving 

(2.1) dX t = fi(X t )dt + a(X t )dB t , 

where fj, and a > are continuous functions satisfying (2.4) and (2.5) be- 
low, and B = (Bt)t>o is a standard Brownian motion. By we denote the 
probability measure under which the process X starts at x > 0. Recalling 
that the scale function of X is given by 

(2.2) L(x)= f expf- [ V - - dz) dy, 



and the speed measure of X is given by 

(2 ' 3) m{dx)= (ay2)t)L'(xY 
we assume that the following conditions are satisfied: 

(2.4) L(0+) = -oo and L(oo-) = 0, 

(2.5) / L(dy) = oo, / m(dy) < oo and / \L(y)\m(dy) < oo. 
Jo+ Jo+ Jo+ 

From (2.4) we read that X is a transient diffusion process in the sense that 
Xt — > oo P^-a.s. as t — > oo, and from (2.5) we read that is an entrance 
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boundary point for X in the sense that the process X could start at but 
will never return to it (implying also that X will never visit after starting 
at x > 0). 

2. The main example we have in mind is the d-dimensional Bessel process 
X solving 

(2.6) dX t = ^-dt + dB u 

where d > 2. Recalling that the scale function is determined up to an affine 
transformation we can choose the scale function (2.2) and hence the speed 
measure (2.3) to read 

(2.7) L(x) = -^ 

2 

(2.8) m(dx) = - — -x^ 1 dx 

for x > 0. It is well known that when d £ {3, 4, . . .} one can realize X as the 
radial part of d-dimensional standard Brownian motion. Similar interpreta- 
tions of (2.6) are also valid when d = 1 (with an addition of the local time 
at zero) and d = 2 but X is not transient in these cases (but recurrent), 
and hence the problem considered below will have a trivial solution. Other 
examples of (2.1) are obtained by composing Bessel processes solving (2.6) 
with strictly decreasing and smooth functions. This is of interest in finan- 
cial applications and will be discussed below. There are also many other 
examples of transient diffusion processes solving (2.1) that are not related 
to Bessel processes. 

3. To formulate the problem to be studied below consider the diffusion 
process X solving (2.1), and introduce its running minimum process / = 
(It)t>0 by setting 

(2.9) I t = inf X s 

0<s<t 

for t > 0. Due to the facts that X is transient (converging to +oo) and 
is an entrance boundary point for X, we see that the ultimate infimum 
/oo = infj>o Xt is attained at some random time 9 in the sense that 

(2.10) X e = Ioo 

with P^-probability one for x > given and fixed (the case x = being triv- 
ial and therefore excluded). It is well known that 9 is unique up to a set 
of P x -probability zero (cf. [42], Theorem 2.4). The random time 9 is clearly 
unknown at any given time and cannot be detected through sequential ob- 
servations of the sample path 1 1— >■ Xt for t > 0. In many applied situations 
of this kind, we want to devise sequential strategies which will enable us to 
come as "close" as possible to 9. Most notably, the main example we have 
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in mind is the problem of optimal trading in the presence of bubbles to be 
addressed below. In mathematical terms this amounts to finding a stopping 
time of X that is as "close" as possible to 9. A first step toward this goal is 
provided by the following lemma. We recall that stopping times of X refer 
to stopping times with respect to the natural filtration of X that is defined 
by = a(X s \0 < s < t) for t > 0. 

Lemma 1. We have 

(2.11) |0-t|=0+ I {21(9<t)-l)dt 

Jo 

for all stopping (random) times r of X . 

Proof. The identity is well known (see, e.g., [35], page 450) and can be 
derived by noting that 

/•O rT 

\0- T \ = (9-t) + + (t-9) + = / I(r<t)dt + I(9<t)dt 

Jo Jo 

= [ (1- I(T>t))dt+ [ I(9<t)dt 

Jo Jo 

rT rT 

(2.12) =9- j I(6>t)dt+ I{6<t)dt 

Jo Jo 

rT rr 

= 0-1 {I- 1(9 <t))dt+ j l(0<t)dt 
Jo Jo 

= 9+1 (21(9<t)-l)dt 
Jo 

for all stopping (random) times r of X as claimed. □ 

4. Taking E x on both sides in (2.11) yields a nontrivial measure of error 
(from r to 9) as long as E x 9 < oo for x > given and fixed. The latter con- 
dition, however, may not always be fulfilled. For example, when X is a tran- 
sient Bessel process of dimension d > 2 it is known (see [38], Lemma 1) that 
P x (9 > t) ~ t~( d / 2 - 1 ) as t -)■ oo. Hence we see that E x 9 = P x (9 >t)dt< 
oo if and only if d/2 — 1 > 1 or equivalently d > 4. It is clear from (2.11), 
however, that the pointwise minimization of the Euclidean distance on the 
left-hand side is equivalent to the pointwise minimization of the integral 
on the right-hand side. To preserve the generality we therefore "normalize" 
\0 — t\ on the left-hand side by subtracting from it. After taking on 
both sides of the resulting identity, we obtain 

(2.13) E x (\9 -t\-9) = E x [ T (2I(9 <t)-l) dt 

Jo 
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for all stopping times r of X (for which the right-hand side is well defined). 
The optimal prediction problem therefore becomes 

(2.14) V(x) = ME x (\9-r\-9), 

T 

where the infimum is taken over all stopping times r of X (with finite mean) 
and x > is given and fixed. Note that the problem (2.14) is equivalent to 
the problem of minimizing E x \9 — r| over all stopping times r of X (with 
finite mean) whenever E x 6 < oo. To tackle the problem (2.14) we first focus 
on the right-hand side in (2.13) above. 

Lemma 2. We have 

(2.15) E x f\2I{6 <t)-l)dt = E x jT (l - 2^y) dt 

for all stopping times r of X ( with finite mean ) and all x > . 

Proof. Using a well-known argument (see, e.g., [35], page 450), we find 
that 

E x / (21(9 < t) - 1) dt = E x / (21(9 <t)- l)I(t < r) dt 
Jo Jo 

/•oo 

= / E X (E X [(2I(9 <t)- l)I(t < r)| Jf ]) dt 
Jo 

roo 

(2.16) = / E x (I(t<r)E x [(2I(9<t)-l)\T t x ])dt 

Jo 

= E X [ T (2P x (9<t\T t x )-l)dt 
Jo 

= E X [ T (I -2P x (9>t\T t x ))dt 
Jo 

for any stopping time r of X (with finite mean) and any x > given and 
fixed. Setting /* = inf s >^ X s and recalling that It = info< s <t X s , we find by 
the Markov property that 

P x (9 > t\T t x ) = P x (I l < I t \T x ) = P x (7* < i\T x )\ l=It 

(2.17) 

= Px(Ioo °9 t < i\JF x )\ i=h = P Xt (Ioo < i)\i=i t 

for t > 0. To compute the latter probability we recall that M := L(X) is a 
continuous local martingale and note that 1^ < i if and only if L(/oo) < L(i) 
where L(/oo) = ^U>o L(X t ) = inft^o-Mj. This shows that the set {1^ < i} 
coincides with the set {mi t >oM t < L(i)} which in turn can be expressed 
as {a < oo} where a = inf{i > 0\Mt < L(i)}. Taking x > i we see that the 
continuous local martingale M" = (M ffAt ) t >o is bounded above by and 
bounded below by L(i) with Mq = L(x) under P x . It follows therefore that 
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M a is a uniformly integrable martingale and hence by the optional sampling 
theorem we find that 

L(x) = E x M a = E x [L(i)I(a < oo)] + E^M^I^a = oo)] 

(2.18) 

= L(i)P x (a <oo) 

upon using that := lim^ooMt = P x -a.s. on {a = oo}. Combining 

(2.18) with the previous conclusions we obtain 

(2.19) P 3 ,(/oo<i) = P*^<oo) = ; ^r 

L(i) 

for i < x in (0,oo). From (2.17) and (2.19) we see that 



(2.20) P x {e>t\F t 



■x, _ L{X t ) 



for all x > and t > (for the underlying three-dimensional law, see [5], 
Theorem A). Inserting this expression back into (2.16) we obtain (2.15) and 
the proof is complete. □ 

5. From (2.13) and (2.15) we see that the problem (2.14) is equivalent to 

(2.21) V(x) = inf E x . £ (l - 2^^) dt, 

where the infimum is taken over all stopping times r of X (with finite mean) 
and x > is given and fixed. Passing from the initial diffusion process X to 
the scaled diffusion process L{X) we see that there is no loss of generality 
in assuming that \i = in (2.1) or equivalently that L(x) = x for x > (with 
L(Xt) —> P^-a.s. as t — >■ oo). Note that the time of the ultimate minimum 
9 is the same for both X and L(X) since L is strictly increasing. Note also 
that r is a stopping time of X if and only if r is a stopping time of L(X). 
To keep the track of the general formulas throughout we will continue with 
considering the general case (when \i is not necessarily zero and L is not 
necessarily the identity function). This problem will be tackled in the next 
section below. 

6. For future reference we recall that the infinitesimal generator of X 
equals 

for x > 0. Throughout we denote T a = m£{t > 0\X t = a} and set r a ^ = T a ATb 
for a < b in (0,oo). It is well known that 

(2.23) P l( v = o) = ^M ^ P.(I,, = t )- tW - tW 



Ub)-L(a) " •'* ' L(b)-L(a) 
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for a < x < b in (0, oo). The Green function of X is given by 
(2.24) 



r ( , {L{b)-L(y)){L{x)-L{a)) 

G a ,b(x,y) = pTTT rr~\ \ia<x<y<b 

L{b) — L{a) 



(L(b) - L(x))(L(y) - L{a)) 
= L(b)-L(a) if«<y<x<6. 

If / : (0, oo) — > M is a measurable function, then it is well known that 
(2.25) E x [ Ta ' b f(X t )dt= [ f(y)G a)b (x,y)m(dy) 

JO J a 

for a < x < b in (0,oo). This identity holds in the sense that if one of the 
integrals exists, so does the other one, and they are equal. 

3. Optimal stopping problem. It was shown in the previous section that 
the optimal prediction problem (2.14) is equivalent to the optimal stopping 
problem (2.21). The purpose of this section is to present the solution to the 
latter problem. Using the fact that the two problems are equivalent this also 
leads to the solution of the former problem. 

In the setting of (2.1)-(2.5) consider the optimal stopping problem (2.21). 
This problem is two-dimensional and the underlying Markov process equals 
(I, X). Setting l\ = i A info< s <t X s for t > enables (I, X) to start at (i, x) 
under P x for i < x in (0,oo), and we will denote the resulting probabil- 
ity measure on the canonical space by Pj jX . Thus under Pj jX the canonical 
process (I,X) starts at (i,x). The problem (2.21) then extends as follows: 

(3.1) V{i,x) = ME i}X [{^- 2 ^P) dt 

for i < x in (0,oo) where the infimum is taken over all stopping times r of 
X (with finite mean). In addition to a and L from (2.1) and (2.2) above, let 
us set 

(3.2) c(i , x ) = i_ 2 ^M 

L(i) 

for i < x in (0,oo). The main result of this section may then be stated as 
follows. 

Theorem 3. The optimal stopping time in problem (3.1) is given by 

(3.3) r* = inf{i > 0|X t >/*(it)}, 

where the optimal boundary /* can be characterized as the minimal solution 
to 

( ~ A) fl( -, = o-\f(i))L'{f(i)) /■/« 4(i,y)[L(y)-L(i)] 

1 • > 1 U c(i, f(i))[L(f(i)) — L(i)] J t a\y)L'{y) V 
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staying strictly above the curve h(i) = L~ l (L{i)/2) for i>0 (in the sense 
that if the minimal solution does not exist, then there is no optimal stopping 
time). The value function is given by 



for i < x < and V(i,x) =0 for x > f*(i) with i>0. 

Proof. 1. It is evident from the integrand in (3.1) that the excursions 
of X away from the running minimum / play a key role in the analysis of 
the problem. In particular, recalling definition (3.2), we see from (3.1) that 
the process (/, X) can never be optimally stopped in the set Co := {(i,x) £ 
S\c(i,x) < 0} where we let S = {(i,x) £ (0,oo) x (0,oo)|i < x} denote the 
state space of the process (I,X). Indeed, if (i,x) £ Co is given and fixed, 
then the first exit time of (I,X) from a sufficiently small ball with the 
centre at (i,x) (on which c is strictly negative) will produce a value strictly 
smaller than (the value corresponding to stopping at once). Defining 



for i > we see that c(i,x) < for x < h(i) and c(i,x) > for x > h(i) 
whenever i < x in (0, oo) are given and fixed. Note that the mapping i \— > h(i) 
is increasing and continuous as well as that h{i) > i for i > with /i(0+) = 
and /i(+oo) = +oo. This shows that Co = {(«, x) G S\i < x < h(i)}. Note in 
particular that Co contains the diagonal {(i,x) £ S\i = x] in the state space. 

2. Before we formalize further conclusions, let us recall that the general 
theory of optimal stopping for Markov processes (see [35], Chapter 1) im- 
plies that the continuation set in the problem (3.1) equals C = {(i,x) £ 
S\V(i,x) < 0}, and the stopping set equals D = {(i,x) £ S\V(i,x) = 0}. It 
means that the first entry time of (I,X) into D is optimal in problem (3.1) 
whenever well defined. It follows therefore that Co is contained in C, and 
the central question becomes to determine the remainder of the set C. Since 
Xt — > oo P x -&.s. as t — > oo it follows that L(Xt) — > P x -a.s. as t — > oo so that 
the integrand in (3.1) becomes strictly positive eventually, and this reduces 
the incentive to continue (given also that the "favorable" set Co becomes 
more and more distant). This indicates that there should exist a point f{i) at 
or above which the process X should be optimally stopped under Pi tX where 
i < x in (0, oo) are given and fixed. This yields the following candidate: 



for an optimal stopping time in (3.1) where the function i i— >■ f(i) is to be 
determined. 



(3.5) 




(3.6) 



h(i) = L-\\L{i)) 



(3.7) 



r f = inf{t > 0\X t > f(I t )} 
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3. Free-boundary problem. To compute the value function V and determine 
the optimal function /, we are led to formulate the free-boundary problem 

(3.8) (L x V)(i,x) = — c(i,x) for i < x < f(i), 

(3.9) V-(i,x)\ x= i + = (normal reflection) , 

(3.10) V(i, x)\ x= fn\_ = (instantaneous stopping) , 

(3.11) Vx(h x )\x=f(i)- = (smooth fit) 

for i > 0, where Lx is the infinitesimal generator of X given in (2.22) above. 
For the rationale and further details regarding free-boundary problems of 
this kind, we refer to [35], Section 13, and the references therein (we note, 
in addition, that the condition of normal reflection (3.9) dates back to [18]). 

4. Nonlinear differential equation. To solve the free-boundary problem 
(3.8)-(3.11) consider the stopping time r/ defined in (3.7) and (formally) 
the resulting function 

(3.12) V f (i,x) = E l>x [ ' c(I t ,X t )dt 

Jo 

for i < x < f(i) in (0, oo) given and fixed. Applying the strong Markov prop- 
erty of (I, X) at T iJ{i) = inf{t > Q\X t £ (t,/(t))} and using (2.23)-(2.25), we 
find that 

L(f(i))-L(x) 

(3.13) V f (i, x) = V f (i, i) L y y _ ^ ' + c{i,y)G Lm {x,y)m{dy). 

It follows from (3.13) that 

„,. L(f(i))-L(i) jr . 

v ^ i)= ^M^w) Vfihx) 

" L(m)-L(x) I <hy)G i>m (x,y)m(dy). 

Using (3.10) and (3.11) we find after dividing and multiplying with x — f{i) 
that 

V f (i,x) 1 dVf 

(3.15) hm T7TTV\ FTT = ~T77TTy\ 

xt/(i) !(/(»)) " L{x) L'(f(i)) dx 



= 0. 

x=m- 



Moreover, it is easily seen by (2.24) that 
(3.16) 



% W§TW) f } c{h y)GlJ{l) ( "' y)m{dy) 



rf(i) 

= c(i,y)[L(y) - L(i)]m(dy). 

Ji 
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Combining (3.14)-(3.16) we see that 

rfd) 

(3.17) Vf (i,i) = - c(i,y)[L(y)-L(i)]m(dy). 



Inserting this back into (3.13) and using (2.24) and (2.25), we conclude that 

rfii) 

(3.18) V f (i,x) = - c(i,y)[L(y) - L(x)]m(dy) 



for i<x< f{i) in (0,oo). Finally, using (3.9) we find that 

a 2 (f(i))L'(f(i)) 



(3.19) 



2c(i,f(i))[L(f(i))-L(i)] 

'/(*) 

c 'i(hy)l L (y) - L{i)]m{dy) 



for i > 0. Recalling that Co is contained in C, we see that there is no re- 
striction to assume that each candidate function / solving (3.19) satisfies 
f(i) > h(i) for all i > 0. In addition we will also show below that all points 
(i,h(i)) belong to C for i > so that (at least in principle) there would be 
no restriction to assume that each candidate function / solving (3.19) also 
satisfies f(i) > h(i) for all i > 0. These candidate functions will be referred 
to as admissible. We will also see below, however, that solutions to (3.19) 
"starting" at h play a crucial role in finding/describing the solution. 

Summarizing the preceding considerations, we can conclude that to each 
candidate function / solving (3.19), there corresponds the function (3.18) 
solving the free-boundary problem (3.8)-(3.11) as is easily verified by direct 
calculation. Note, however, that this function does not necessarily admit 
the stochastic representation (3.12) (even though it was formally derived 
from this representation). The central question then becomes how to select 
the optimal boundary / among all admissible candidates solving (3.19). To 
answer this question we will invoke the subharmonic characterization of the 
value function (see [35], Chapter 1) for the three-dimensional Markov process 
(I, X, A) where A t = f c(I s ,X s ) ds for t > 0. Fuller details of this argument 
will become clearer as we progress below. It should be noted that among 
all admissible candidate functions solving (3.19) only the optimal boundary 
will have the power of securing the stochastic representation (3.12) for the 
corresponding function (3.18). This is a subtle point showing the full power 
of the method (as well as disclosing limitations of the optimal stopping 
problem itself). 

5. The minimal solution. Motivated by the previous question we note 
from (3.18) that / 1— > Vf is decreasing over admissible solutions to (3.19). 
This suggests to select the candidate function among admissible solutions 
to (3.19) that is as far as possible from h. The subharmonic characterization 
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of the value function suggests to proceed in the opposite direction, and this 
is the lead that we will follow in the sequel. 

To address the existence and uniqueness of solutions to (3.19), denote the 
right-hand side of (3.19) by &(i,f(i)). From the general theory of nonlinear 
differential equations, we know that if the direction field (i,f) i-> /) is 
(locally) continuous and (locally) Lipschitz in the second variable, then the 
equation (3.19) admits a (locally) unique solution. For instance, this will be 
the case if, along a (local) continuity of (i, f) *— > $(i, /), we also have a (local) 
continuity of (i, f) t— > &'Ai, /)• In particular, we see from the structure of $ 
that equation (3.19) admits a (locally) unique solution whenever x i— > a 2 (x) 
is (locally) continuously differentiable. It is important to realize that the 
preceding arguments apply only away from h since each point (i,h(i)) is 
a singularity point of equation (3.19) in the sense that f'(i+) = oo when 
f(i+) = h(i) due to c(i, h(i)) = for i > 0. In this case it is also important to 
note that the preceding arguments can be applied to the equivalent equation 
for the inverse of i t— > f(i) since this singularity gets removed (the derivative 
of the inverse being zero). 

To construct the minimal solution to (3.19) staying strictly above h, we 
can proceed as follows (see Figure 1). For any i n > such that i n \, as 



i i = x 








Fig. 1. Solutions f n and /» to the differential equation (3.4) from Theorem 3. The 
optimal stopping boundary /» is the minimal solution staying strictly above the curve h. 
This is a genuine drawing corresponding to the golden ratio rule of Corollary 5 when X is 
the radial part of three-dimensional Brownian motion and the optimal stopping boundary 
f„ is linear. 
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n — > oo, let i i— >■ f n {i) denote the solution to (3.19) on (i n ,oo) such that 
fn(in+) = h(i n ). Note that i i— > f n (i) is singular at i n and that passing to the 
equivalent equation for the inverse of i i— >■ f n (i), this singularity gets removed 
as explained above. (Note that the solution to the equivalent equation for 
the inverse can be continued below h(i n ) as well until hitting the diagonal 
at some strictly positive point at which the derivative is — oo. This yields 
another solution to (3.19) staying below f n and providing its "physical" 
link to the diagonal. We will not make use of this part of the solution in 
the sequel.) Note that the right-hand side of equation (3.19) is positive 
for f(i) > h(i) so that i h-> f n (i) is strictly increasing on [%, oo). By the 
uniqueness of the solution we know that the two curves i i— > f n (i) and i i— > 
f m (i) cannot intersect for n^m, and hence we see that (/ n )n>i is increasing. 
It follows therefore that /* := lim^-j.oo f n exists on (0, oo). Passing to an 
integral equation equivalent to (3.19) it is easily verified that i h-> f*(i) solves 
(3.19) wherever finite. This /* represents the minimal solution to (3.19) 
staying strictly above the curve h on (0,oo). We will first consider the case 
when /* is finite valued on (0, oo). 

6. Stochastic representation. We show that the function (3.18) associated 
with the minimal solution /* admits the stochastic representation (3.12). For 
this, let i i — y f n {i) be the solution to (3.19) on (i n ,oo) such that f n {i n +) = 
h(i n ) for i n > with i n 1 as n — > oo. Consider the function (z, x) h-» Vf n (i, x) 
defined by (3.18) for i < x < f n (i) and i > i n with n > 1 given and fixed. 
Recall that Vf n solves the free-boundary problem (3.8)-(3.11) for i > i n . 
Consider the stopping time r n := Tj n A rj n where Tj n = inf{i > 0\X t = i n } 
and Tj n = inf{t > 0\X t > f n (It)}. Applying Ito's formula and using (3.8), we 
find that 



V fn (I Tnt X Tn ) = V fn (i,x) + 



I 

Jo 



Tn dV f , 



71 



(It,X t )dI t 



di 




(3.20) 
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where we also use (3.9) to conclude that the integral with respect to dlt 
is equal to zero and M t = /q At " a(X s )(dVf n /dx)(I s , X s ) dB s is a continuous 
local martingale for t > 0. 

Since the process (I,X) remains in the compact set {(j,y) £ S\i n <j < 
V < fn{i)} up to time r n under Pi >x , and both o and dVf n /dx are continuous 
(and thus bounded) on this set, we see that M is a uniformly integrable 
martingale, and hence by the optional sampling theorem we have Ei >x M Tn = 
0. Taking Ei )X on both sides of (3.20), we therefore obtain 



(3.21) 



V fn (i, x) = E i>x V fn (I Tn , X Tn ) + j " c(I t ,X t )dt 

Jo 

f-T n 

= V fn (i n ,i n ) Pi, x (n n < TfJ + E i>x / c(I t ,X t 

Jo 



dt 



since (I Tn ,X Tn ) = (i n , i n ) on {r in < r fn } and V fn (I Tn ,X Tn )=0 on {r /n < r in } . 
Using that |c| < 1 we find by (3.17), (2.23) and (3.6) that 

\Vf n (i n ,in)\ Pi,x( T in < T /J 
/•/(in) 

< / \c(i,y)\\L(y) - L(i n )\m(dy)P iiX (T in <T fn ) 
(3.22) Jin 

<\L(h(i n ))-L(i n )\ f hM m(dy) L{fS)) ~ L{x) 



I n 



L{U(i)) - L{i n ) 



l |r , s . L(U{i))-L{x) f h M 

as n — )• oo since L(i n ) — >■ — oo and /i(i„.) — > 0, so that jf 1 m(dy) — > 0, due 
to (2.5) above. Hence letting n — > oo in (3.21) and using that Vf n — > Vf, by 
the monotone convergence theorem, as well as that r n f since / n f /* and 
i n 4, 0, we find noting that Ei^Tf^ < 00 and using the dominated convergence 
theorem that 

(3.23) V f , (i, x) = E i>x [ U c(I t , X t ) dt 

Jo 

for all % < x in (0, 00) as claimed. 

7. Nonpositivity. We show that for every solution / to (3.19) such that 
/ > /* on (0,oo) and the function Vf defined by (3.18) above, we have 

(3.24) V f (i,x)<0 

for all i < x in (0,oo). Clearly, since c(i,y) > for y > h{i) in (3.18), it is 
enough to prove (3.24) for and i <x < h(i) with i > 0. For this, consider 
the stopping time = inf{i > Q\X t > h(It)} and note that r/, = + ry, o 
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9 Th . Hence by the strong Markov property of (I,X) applied at Th we find 
using (3.23) that 



Th 



c(I t ,X t )dt + Ei tX 



Th+T U oU Th 



(3.25) 



Ej x 



c(I t ,X t )dt + E itX 



T f* O0 T h 







c(I t ,X t )dt 
c{I t+Th ,X t+Th )dt 



Th 



c(I t ,X t )dt + E ijX E itX 



c(I t ,X t )dtoO 



Th 



Th 



c(I t ,X t ) dt + Ei iX E/ >x 



c(I t ,X t )dt 



= E i>x I h c(I t ,X t )dt + E i>x V ft (I Th ,X Th )< 0, 
Jo 

where the final inequality follows from the facts that c(It,Xt) < for all 
t e [OjTh] and Vf t (I Th , X Th ) < due to X Th = h(I Th ) upon recalling (3.18) as 
already indicated above. This completes the proof of (3.24). 

8. Optimality of the minimal solution. We will begin by disclosing the 
subharmonic characterization of the value function (3.1) in terms of the 
solutions to (3.19) staying strictly above h. For this, let i i— > f(i) be any 
solution to (3.19) satisfying f{i) > h(i) for all i > 0. Consider the function 
(i, x) i — y Vf(i, x) defined by (3.18) for i < x < f(i) in (0, oo) and set Vf(i, x) = 
for x > f(i) in (0, oo). Let i < x in (0, oo) be given and fixed. Due to the 
"double-deck" structure of Vf, we can apply the change-of- variable formula 
from [32] that in view of (3.11) reduces to standard Ito's formula and gives 

f* dV f 

Vf(I t ,X t ) = V f (i,x) + -J-(I s ,X s )dI s 

/"* dVt 1 f* d 2 V f 

+ J -d {Is > Xs)dXs + 2 I -^-( I s,X s )d(X,X) s 

= V f (i,x)+ f h x (V f )(I s ,X s )ds 
Jo 



(3.26) 



+ 



o 



' 1 a(X s )^(I s ,X s )dB s 



where we also use (3.9) to conclude that the integral with respect to dl s is 
equal to zero. The process M = (M t )t>o defined by 



(3.27) 



M t 



dV f 

a{X s )-^(I s ,X s )dB s 
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is a continuous local martingale. Introducing the increasing process P = 
(Pt)t>0 by setting 

(3.28) P t = f c{I s ,X s )I(X s > /(/,)) da 

Jo 

and using the fact that the set of all s for which X s equals f(I s ) is of 
Lebesque measure zero, we see by (3.8) that (3.26) can be rewritten as 
follows: 

(3.29) Vf{It,X t )+ I c{I s ,X s )ds = V f (i,x)+M t + P t . 

Jo 

Prom this representation we see that the process Vf(It,X t ) + J * c(J s , X s ) ds 
is a local submartingale for t > 0. 

Let r be any stopping time of X (with finite mean) . Choose a localization 
sequence (<7 n )n>i of bounded stopping times for M. Then by (3.24) and 
(3.29) we can conclude using the optional sampling theorem that 

rrAa n rTAcr„ 

E ijX c{I t ,X t )dt> E i>x V f (I TAan ,X TAan )+ / c(I t ,X t )dt 

Jo l Jo 

(3.30) 

> V f (i,x) + E i)X M rAan = V f (i,x). 

Letting n — > oo and using the dominated convergence theorem (upon recall- 
ing that |c| < 1 as already used above) we find that 

(3.31) E hX [ c(I t ,X t )dt>V f (i,x). 

Jo 

Taking first the infimum over all r, and then the supremum over all /, we 
conclude that 

(3.32) V{i,x)>supV f {i,x) = V ft (i,x), 

f 

upon recalling that / 1— > Vf is decreasing over / > /* so that the supremum 
is attained at /*. Combining (3.32) with (3.23) we see that (3.3) and (3.5) 
hold as claimed. 

Note that (3.30) implies that the function (i,x) \— > Vj(i,x) + a is subhar- 
monic for the Markov process (I, X, A) where A t = £ c(I s ,X s ) ds for t > 0. 
Recalling that / 1->- Vf is decreasing over / > /*, and that Vf(i,x) < for all 
i < x in (0, oo) by (3.24) above, we see that selecting the minimal solution 
/* staying strictly above h is equivalent to invoking the subharmonic char- 
acterization of the value function (according to which the value function is 
the largest subharmonic function lying below the loss function). For more 
details on the latter characterization in a general setting we refer to [35], 
Chapter 1. It is also useful to know that the subharmonic characterization of 
the value function represents the dual problem to the primal problem (3.1) 
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(for more details on the meaning of this claim including connections to the 
Legendre transform see [33]). 

Consider finally the case when /* is not finite valued on (0,oo). Since 
i i — y f*{i) is increasing we see that there is i* > such that /*(i) < oo for all 
i £ (0, i*) when > and /*(i) = oo for all i > i* with i ^ when i* = 0. If 
i* > 0, then the proof above can be applied in exactly the same way to show 
that (3.3) and (3.5) hold as claimed under Pj jX for all i < x in (0,oo) with 
i <i*. If i > i* with i ^ when = 0, then the same proof shows that (3.5) 
still holds with oo in place of /*(«), however, the stopping time (3.3) can 
no longer be optimal in (3.1). This is easily seen by noting that the value 
in (3.5) is nonpositive (it could also be — oo) for any x > h(i) for instance, 
while the Pj iX -probability for X hitting i before drifting away to oo is strictly 
smaller than 1 so that the Pj^-expectation over this set in (3.1) equals oo 
(since the integrand tends to 1 as t tends to oo) showing that the stopping 
time (3.3) cannot be optimal. The proof above shows that the optimality 
of (3.5) in this case is obtained through r n = Tj n A ry n which play the role of 
approximate stopping times (obtained by passing to the limit when n tends 
to oo in (3.21) above). This completes the proof of the theorem. □ 

4. The golden ratio rule. In this section we show that the minimal so- 
lution to (3.4) admits a simple closed-form expression when X is a tran- 
sient Bessel process (Theorem 4). In the case when X is the radial part 
of three-dimensional Brownian motion this leads to the golden ratio rule 
(Corollary 5). We also show that X stopped according to the golden ratio 
rule has what we refer to as the golden ratio distribution (Corollary 8). 

In the setting of (2.6)-(2.8) consider the optimal prediction problem 
(2.14). Recall that this problem is equivalent to the optimal stopping prob- 
lem (2.21) which further extends as (3.1). The main result of this section 
can now be stated as follows. 

Theorem 4. If X is the d- dimensional Bessel process solving (2.6) with 
d>2, then the optimal stopping time in (2.14) is given by 

(4.1) r* = inf{t > 0\X t > XI t }, 
where A is the unique solution to 

(4.2) x d -(l + d)\ 2 + -^\ A - d -^f^ = Q ifd^A, 

A-d 4-d 

(4.3) A 4 -5A 2 + 41ogA + 4 = ifd = A 
belonging to (2 1 ^ d ~ 2 ),oo). The value function (3.1) is given explicitly by 
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(4.4) 



x 

~ ~d 



2X 



4-d 



d-4 



d-4 



ijd + 4 



-I 

x 



2? log I ^ 



i/d = 4 



/or i < x < Xi and V(i, x) = /or x> Xi with i > 0. 

Proof. By the result of Theorem 3 we know that the optimal stopping 
time r* is given by (3.3) above where the optimal boundary /* can be char- 
acterized as the minimal solution to (3.4) staying strictly above the curve 
h(i) = L~ 1 (L(i)/2) for i > 0. Using (2.7) and (3.2) it can be verified that 
(3.4) reads as follows: 

d-2 / £ i -w d-4 



f'ii) 
(4.5) 



d-2 



4-dV 



(m 



(4-d) 



d-2 



+ (d - 2) 



d-2 



if d^4 



2(/(z)A)[(/(i)A) 2 -21og(/(z)/i)-l] 



if d = 4 



((/(<)/i) 2 -l)((/«/<) 2 -2) 

and /i(i) = 2 1 /( d ~ 2 )i for i > 0. Hence it is enough to show that f*(i) = Xi is 
the minimal solution to (4.5) staying strictly above the curve h(i) = 2 1 ^ d ~ 2 ^i 
for i > 0. 

To show that /* is a solution to (4.5) staying strictly above h, insert 
f(i) = Xi into (4.5) with A > to be determined. Multiplying both sides of 
the resulting identity by X A ~ d (to be able to derive the factorization (4.7) 
below) it is easy to see that this yields the equation F(X) = where we set 

(d-2) 2 



F (X) = X a -(l + d)X' + 



-X 



4-d 



(4.6) 



4-d 



4-d 

A 4 -5A 2 +41ogA + 4 ifd = 4 
for A > 0. After some algebraic manipulations we find that 

2 



if d^4 



(4.7) 



F'(X)=dX 3 - d { X d ~ 2 



d 



(A 



d-2 



2) 



for A > and d > 2. Hence we see that the equation F'(X) = has two roots 
A = (2/d) 1 /( d ~ 2 ) and Ai = 2 1 /^~ 2 ) where < A < 1< Ai < oo. It is easy to 
check that F"(Xq) < and F"(X\) > showing that F has a local maximum 
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at Ao and F has a local minimum at Ax- Noting that F(0+) < 0, -F(l) = 
and F(oo— ) = oo this shows that (i) F is strictly increasing on (0, Ao) 
with F(0+) < and F(Aq) > 0; (ii) F is strictly decreasing on (Ao, Ai) with 
F(l) = and F(\i) < 0; and (hi) F is strictly increasing on (Ax,oo) with 
F(oo— ) = oo. It follows therefore that the equation F{X) = has exactly 
three roots Aq < 1 < A^ where Aq G (0, Ao) and A^ G (Ai,oo). Setting A = A^ 
this shows that = Ai is a solution to (4.5) staying strictly above the 
curve h(i) = 2 1 /( d ~ 2 H for i > as claimed. 

To show that /* is the minimal solution satisfying this property, set = 
f(i)/i and note that (4.5) can then be rewritten as follows: 

(4.8) IK (l) = ~ K Z-d^ K d-2^ _ l)( K d-2(j) _ 2 ) 

for i > 0. Since F{n{i)) < for G (2 1 /( d_2 ),A) we see from (4.8) that 
z I—)- «(i) is increasing for k(i) G (2 1 /( rf_2 ^, A). Noting that (4.8) implies that 

(4.9 - / dn= — = log — 

it follows therefore that the integrand on the left-hand side is bounded by 
a constant (not dependent on i) as long as «(t) G (2 1 /( d ~ 2 ),A) for i G (0,i ) 
with any iq > given and fixed. Letting then i 1 in (4.9) we see that 
the left-hand side remains bounded while the right-hand side tends to oo 
leading to a contradiction. Noting that G {2 1 ^ <yd ~ 2 \\) if and only if 
f(i) G (h(i), f*(i)), we can therefore conclude that there is no solution / to 
(4.5) satisfying f(i) G (h(i), f*(i)) for i > 0. Thus /* is the minimal solution 
to (4.5) staying strictly above h and the proof is complete. □ 

Corollary 5 (The golden ratio rule). If X is the radial part of three- 
dimensional Brownian motion, then the optimal stopping time in (2.14) is 
given by 

Xt-h 



(4.10) r* = inNi>0 



It 

where ip = (1 + v / 5)/2 = 1.61 . . . is the golden ratio (see Figure 2). 

Proof. In this case d = 3 and equation (4.2) reads 

(4.11) A 3 - 4A 2 + 4A - 1 = (A - 1)(A 2 - 3A + 1) = 

for A > 0. Solving the latter quadratic equation and choosing the root strictly 
greater than 1, we find that A = 1 + <p where ip = (1 + v5)/2 = 1.61 ... is the 
golden ratio. The optimality of (4.10) then follows from (4.1) and the proof 
is complete. □ 
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1 t^{l + <fi)I t 

1 




t 



Fig. 2. The golden ratio rule for the radial part X of three-dimensional Brownian motion. 

Returning to the result of Theorem 3 above we now determine the law of 
the transient diffusion process X stopped at the optimal stopping time r* 
(for related results on the Skorokhod embedding problem see [36], pages 269- 
277, and the references therein). 

Proposition 6. In the setting of Theorem 3 we have 

(4.12) P.(Xr.<v) = exp(- f — f* z) ) 
for < y < f*(x) with x > 0. 

Proof. Note that 

n=mf{t>0\X t >/*(/*)} 

(4.13) = m£{t > 0\L(X t ) >(io/,o L-^mit))} 
= mf{t>0\X^>f^I t L )}, 

where we set = L(X t ), f£ = Lo/ t oL _1 and J/ 1 = L(I t ) = inf < s <t L(X S ) = 
info< s <tXf for t > 0. Let x > be given and fixed. For j < L(x) set G(j) = 
flocditydk where g: (— oo,L(x)] — > R is a continuously differentiable func- 
tion with bounded support. Using the fact that dl^ = when X^ ^ 1^ it is 
easily verified by Ito's formula that the process M L = (M^)t>o defined by 

(4.14) M t L = G(I t L ) + (X t L - I t L )G'(I t L ) 

is a continuous local martingale. Moreover, since G' = g is continuous and 
has bounded support, we see that M L is bounded and therefore uniformly 
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integrable. By the optional sampling theorem we thus find that 
G L (x) = E X M£ = E X M^ = E X G(I^) + E X [(X^ - 

= [ L(X) G(j) dF(j) + E x [(f^J - 

J — oo 

= G(j)F(j)\ L _% - [ L(X) F(j) dGU) 

(4.15) 

LX \fHj)-j)G'(j)dF(j) 

-oo 

r L(x) 

G L (x)- / F(j)g(j)dj 

J — oo 

L{X \fHj)-j)g(j)dF{j), 

-oo 

where we set G L (x) = G{L{x)) and F denotes the distribution function of 
1% under P x . Since (4.15) holds for all functions g of this kind, it follows 
that 

(4 ' 16> F ' ( ^7FIJF7 

for j < L(x) with F(L(x)) = 1. Solving (4.16) under this boundary condition 
we find that 

/ rL{x) \ 

for j < L(x). Recalling that = L o o L^ 1 and substituting k = L(z) it 
follows that 

P*(Ir. < *) = Px(i(IrJ < L(i)) = P X (I^ < L(i)) 

/ r L ( x ) dk \ 

(4 - 18) = Fim M-L 7m^) 

dL(z) 



exp 



L{U(z))-L(z) 



for i < x in (0, oo). Hence we find that 

P X {X T , <y) = P*(/,(/ T .) < y) = P»(/ T . < Z" 1 ^)) 

(4.19) 

■ exp' 



f-\ y) L(f m (z))-L(z) 
for < y < f*(x) with 3; > 0. This completes the proof. □ 



22 



K. GLOVER, H. HULLEY AND G. PESKIR 



Specializing this result to the ci-dimensional Bessel process X of Theo- 
rem 4 we obtain the following consequence. 

Corollary 7. In the setting of Theorem 4 we have 

(d-2)/(l-(l/A) d - 2 ) 

(4.20) Px{X T ,<y) 
for < y < Xx with x > 0. 



Xx) 



Proof. In this case /*(&) = Xi for i > where A is the unique solution to 
either (4.2) when d ^ 4 or (4.3) when d = 4 and L is given by (2.7). Inserting 
these expressions into the right-hand side of (4.12) it is easily verified that 
this yields (4.20). □ 

Specializing this further to the radial part X of three-dimensional Brow- 
nian motion in Corollary 5 we obtain the following conclusion. 

Corollary 8 (The golden ratio distribution). In the setting of Corol- 
lary 5 we have 

(4.21) P x (X Tf <y) = 



y 



[l+<p)x 



for < y < (1 + ip)x with x > 0. 



Proof. In this case d = 3 and A = l + (p so that (d-2)/(l - (1/A) d ~ 2 ) = 
l/(l — l/(l + i^)) = (l + ip)/(p = tp 2 /if = if. Hence we see that (4.20) reduces 
to (4.21). □ 

Note from (4.21) that the density function of X Tt under P x is given by 

for < y < (1 + (p)x with x > and equals zero otherwise. We refer to (4.21) 
and (4.22) as the golden ratio distribution. It is easy to see that 

(4.23) E x X Tt = ipx 

for x > 0. The fact that this number is strictly greater than x (the initial 
point corresponding to stopping at once) is not surprising since X is a sub- 
martingale. It needs to be recalled moreover that the aim of applying the 
golden ratio rule r* is to be as close as possible to the time 9 at which 
the ultimate minimum is attained. We will see in the next section that the 
golden ratio distribution provides insight as to what extent the golden ra- 
tio rule has the power of capturing the ultimate maximum of a strict local 
martingale. 



3D BROWNIAN MOTION AND THE GOLDEN RATIO RULE 



23 



5. Applications in optimal trading. In this section we present some ap- 
plications of the previous results in problems of optimal trading. We also 
outline some remarkable connections between such problems and the prac- 
tice of technical analysis. These applications and connections rest on three 
basic ingredients that we describe first. 

1. Fibonnaci retracement. We begin by explaining a few technical terms 
from the field of applied finance. Technical analysis is a financial term used to 
describe methods and techniques for forecasting the direction of asset prices 
through the study of past market data (primarily prices themselves plus 
the volume of their trade). Support and resistance are concepts in techni- 
cal analysis associated with the expectation that the movement of the asset 
price will tend to cease and reverse its trend of decrease/increase at cer- 
tain predetermined price levels. A support / resistance level is a price level at 
which the price will tend to find support/resistance when moving down/up. 
This means that the price is more likely to bounce off this level rather than 
break through it. One may also think of these levels as turning points of 
the prices. Fibonacci retracement is a method of technical analysis for de- 
termining support and resistance levels. The name comes after its use of 
Fibonacci numbers F n+ \ = F n + F n _\ for n > 1 with Fq = and F\ = 1. 
Fibonacci retracement is based on the idea that after reversing the trend 
at a support /resistance level, the price will retrace a predictable portion 
of the past downward/upward move by advancing in the opposite direction 
until finding a new resistance/support level, after which it will return to 
the initial trend of moving downwards/upwards. Fibonacci retracement is 
created by taking two extreme points on a chart showing the asset price as a 
function of time and dividing the vertical distance between them by the key 
Fibonacci ratios ranging from 0% (start of the retracement) to 100% (end of 
the retracement representing a complete reversal to the original trend). The 
other key Fibonacci ratios are 23.6% (shallow retracement), 38.2% (mod- 
erate retracement) and 61.8% (golden retracement). They are obtained by 
formulas (F n /F n+3 ) x 100 « ip~ 3 x 100, (F n /F n+2 ) x 100 ps ip~ 2 x 100 and 
(F n /F n+ i) x 100 w tp^ 1 x 100, respectively (see the next paragraph). These 
retracement levels serve as alert points for a potential reversal at which 
traders may employ other methods of technical analysis to identify and con- 
firm a reversal. Despite its widespread use in technical analysis of asset 
prices, there appears to be no (rigorous) explanation of any kind as to why 
the Fibonacci ratios should be used to this effect. We will show below that 
the golden ratio rule derived in the previous section offers a rigorous opti- 
mality argument for the choice of the golden retracement (61.8%). To our 
knowledge this is the first time that such an argument has been found/given 
in the literature. 
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2. Golden ratio and Fibonacci numbers. The link between the two is well 
known and is expressed by Binet's formula 

, x ip n -tb n m n -lb n 

(5.1) F, - 



n 



<P~ip y/5 

where <p=(l + Vb) /2 and i/> = (1 - y/l)/2 = 1 - <p = -l/<p. It follows that 

(5.2) lim = 

This fact is used in the description of Fibonacci retracement above. 

3. The CEV model. One of the simplest/tractable models for asset price 
movements that is capable of reproducing the implied volatility smile/frown 
effect and the (inverse) leverage effect (both observed in the empirical data) 
is the Constant Elasticity of Variance (CEV) model in which the (nonneg- 
ative) asset price process Z = (Z t )t>o solves 

(5.3) dZ t = ^iZ t dt + aZl +p dB u 

where \x G M. is the appreciation rate, a > is the volatility coefficient, and 
j3 G R is the elasticity parameter. If (3 = then Z is a geometric Brown- 
ian motion which was initially considered in [28] and [37]. For /3 ^ this 
model was firstly considered in [4] for j3 < and then in [13] for /3 > 0. 
Due to its predictive power and tractability, the CEV model is widely used 
by practitioners in the financial industry, especially for modeling prices of 
equities and commodities. If f3 < then the model embodies the leverage 
effect (commonly observed in equity markets) where the volatility of the 
asset price increases as its price decreases. If @ > then the model embod- 
ies the inverse leverage effect (often observed in commodity markets) where 
the volatility of the asset price increases when its price increases. For ex- 
ample, it is reported in [17] that the elasticity coefficient (3 for Gold on the 
London Bullion Market in the period from 2000 to 2007 was approximately 
0.49. Similar elasticity coefficients have also been observed for other precious 
metals (such as Copper for instance). 

In the remainder of this section we focus on the case when fi = and /3 > 0. 
It is well known (cf. [13]) that Z solving (5.3) is a strict local martingale (a 
local martingale which is not a true martingale) in this case due to the fact 
that 1 1 — y E z (Zt) is strictly decreasing on K + for any z > 0. This also implies 
that Z does not admit an equivalent martingale measure so that the CEV 
model may admit arbitrage opportunities. One way of looking at the models 
of this type is to associate them with asset price bubbles (see [21]). After 
soaring to a finite ultimate maximum (bubble) at a finite time, the asset 
price will tend to zero as time goes to infinity, and the central question for 
a holder of the asset becomes when to sell so as to be as close as possible 



3D BROWNIAN MOTION AND THE GOLDEN RATIO RULE 



25 



to the time at which the ultimate maximum is attained. More precisely, 
introducing the running maximum process S = (St)t>o associated with Z 
by setting 

(5.4) St = sup Z s 

0<s<t 

and recalling that Zt — > as t — > oo, we see that the ultimate supremum 
Soo = sup t>0 Zt is attained at some random time 6 in the sense that 

(5.5) Z g = 

with P^-probability one for z > given and fixed. The optimal selling prob- 
lem addressed above then becomes the optimal prediction problem 

(5.6) V(z) = mfE z (\9-r\-6), 

T 

where the infimum is taken over all stopping times r of Z (with finite mean) 
and z > is given and fixed. We will now show that due to the well-known 
connection between CEV and Bessel processes (dating back to similar trans- 
formations in [4] and [13]) the problem (5.6) can be reduced to the problem 
(2.14) solved above. 

4. The golden ratio rule for the CEV process. For d > 2 given and fixed 
consider the d-dimensional Bessel process X solving (2.6) under P x with x > 
0. Recall that the scale function L is given by (2.7) and set K(x) = —c a L(x) 
for x > with c CT > given and fixed. Then the process Z = K(X) defined 
by 

(5-7) Z t = K(X t ) = ^ 

is on natural scale and Ito's formula shows that Z solves 
(5.8) dZ t = aZl +1/(d - 2) dB t , 

where a = (d — 2) jcj^ 2 ^ and B = — B is a standard Brownian motion. Note 
that equation (5.8) coincides with equation (5.3) for fi = and j3 = l/(d — 2). 
By the uniqueness in law for this equation (among positive solutions) it 
follows that Z = K{X) is a CEV process. From the properties of X it follows 
that after starting at z = K{x) > 0, the process Z stays strictly positive 
(without exploding at a finite time) and Z% — > with P^-probability one 
as t — > oo. This shows that 6 in (5.5) is well defined. Moreover, due to the 
reciprocal relationship (5.7) we see that the time of the ultimate maximum 
6 for Z in (5.5) coincides with the time of the ultimate minimum 8 for X in 
(2.10), and hence the problem (5.6) has the same solution as the problem 
(2.14) (note also that the natural nitrations of Z and X coincide so that r is 
a stopping time of Z if and only if r is a stopping time of X). Since Xt > Xlt 
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Fig. 3. The golden ratio rule for the CEV process Z = 1/X where X is the radial part 
of three-dimensional Brownian motion. Note the presence of a bubble and its relation to 
the qolden ratio. 



if and only if Zt/c a 



X?~ d < A 



2-d T 2-d 



Y~ d St/c a it follows from (4.1) in 



Theorem 4 that the optimal stopping time in (5.6) is given by 
(5.9) r* = inf{* > 0\S t > X d ' 2 Z t }, 

where A is the unique solution to either (4.2) or (4.3) belonging to (2 1 / ( d-2 ) , oo) 



In particular, if d = 3 then we know from (4.11) that A 
reads 



1 + ip so that (5.9) 



(5.10) 



r* = inf<U>0 



St-Zt 



This is the golden ratio rule for the CEV process Z = 1 /X where X is the 
radial part of three-dimensional Brownian motion (see Figure 3). 

To relate the golden ratio rule (5.10) to Fibonacci retracement discussed 
above, let a = S Tt — Z Tt denote the larger quantity and let b = Z Tf denote 
the smaller quantity in the golden ratio rule. To determine the percentage 
of a in a + b we need to calculate 



St* ~ Z Tf 



(5.11) 



a + b 



1 



1 



1 



1 



Multiplying this expression by 100 gives 61.8% and this is exactly the golden 
retracement discussed above. In view of the optimality of (5.10) in (5.6) we 
see that the golden retracement of 61.6% for the CEV process Z = 1/X 
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(starting close to zero) where X is the radial part of three-dimensional Brow- 
nian motion can be seen as a rational support level (in the sense that rational 
investors who aim at selling the asset at the time of the ultimate maximum 
will sell the asset at the time of the golden retracement, and therefore the 
asset price could be expected to raise afterwards). To our knowledge this 
is the first time that such a rational optimality argument for the golden 
retracement has been established. 
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